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Abstract 

A dripline nucleus 22 C is studied in a Borromean three-body model of 20 C+n+n. The valence 
neutrons, interacting via a realistic potential, are constrained to be orthogonal to the occupied 
orbits in 20 C. We obtain ample results supporting that 22 C is an ideal S"-wave two-neutron halo 
nucleus: The ground state is bound by 390-570 keV, the root mean square neutron and proton 
radii are 4.0 and 2.4 fm, and the two neutrons are predominantly in (si/2) 2 orbits. The binding 
mechanism of 22 C is discussed. One- and two-body density distributions elucidate the halo char- 
acter as well as the correlated motion of the neutrons. The reaction cross sections of 22 C+ 12 C 
collisions are predicted. 
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The subshell closure of N=14 and iV=16 is one of the topics discussed intensively in 
the study of neutron-rich nuclei, and the N=14 closure has experimentally been confirmed 
around 22 |l|, |2j, La, \% |5( . This issue is closely related to the competition of 0d 5 / 2 and lsi/2 
neutron orbits. In fact, they play a vital role in determining the ground state structure of 
A=15-20 carbon isotopes. For example, the ground state of 16 C is found to contain the 
(lsi/2) 2 and (Oc/5/2) 2 configurations nearly equally .whereas the last neutron in 19 C is 
in the lsi/2 orbit, forming one-neutron halo structure l a |9j|. The deformations of carbon 
isotopes are discussed to have a strong iV-dependence |l(J,lll|]. An overview of the structure 
of carbon isotopes is given in Ref. |l2j]. 

No information is available to determine whether the iV=14 subshell closure occurs in 
20 C. The systematics of the interaction cross section suggests, however, that the radius 
of 20 C is smaller than that of 19 C ^3|, so it is natural to assume that the ground state 
of 20 C predominantly consists of a (0<i5/2) 6 configuration. If its dominant component were 
(0rf 5/ /2) 4 (lsi/2) 2 , one more neutron could be added to the 0^5/2 orbit to form a particle-stable 
21 C, which is in contradiction to observation. 

In this study we will demonstrate that 22 C is an s-wave two-neutron halo nucleus on the 
basis of the analysis of its structure including the neutron and proton densities. For Z < 8, 
22 C is an only dripline nucleus which the interaction or reaction cross section measurement 
on a 12 C target has not reached yet [l3j . so it is of interest for a future measurement to 
predict the reaction cross section of 22 C. The neutron and proton densities obtained here 
will be useful to estimate the reaction cross section of 22 C on a proton target, which is being 
investigated experimentally [14j. Our model is that 22 C is a three-body system of 20 C+n+n, 
and that 20 C has the (Og^) 6 configuration. 22 C is thus Borromean, just as n Li is. Though 
22 C may be expected to be much like n Li in its halo character, a remarkable difference 
will show up: In n Li both of (Opi/2) 2 and (lsi/2) 2 components contribute to producing its 
halo (liL Ha. [rj, whereas in 22 C only an (si/2) 2 component will be predominant. Another 
difference to be noted is that the 20 C core has zero spin, which will make the content of 
angular momentum coupling in 22 C simpler than that in 11 Li. 

The wave function for 22 C is determined from the following Hamiltonian 

H = T x + T p + U 1 + U 2 + v X2 , (1) 

where the subscripts, A and p, of the kinetic energies stand for the relative distance vectors 
of the three-body system. The two-neutron potential v 12 is taken from the realistic G3RS 
(case 1) potential |l8( which contains central, tensor and spin-orbit forces and reproduces 
the nucleon-nucleon scattering data as well as the deuteron properties. Ui is the n— 20 C 
potential whose form is assumed as 

U = -V f(r) + V x l ■ 8-^-f(r) + V s e~^V s , (2) 

r ar 

where f(r) = [l + exp( r ~ g Rc )]~ 1 with R c =r A^ 3 (A c = 20). The operator V s of the last term 
projects to the s wave of the n— 20 C relative motion, so this term modifies the s-wave 
potential strength. To determine the parameters of U, we take into account the conditions 
that (i) the lsi/2 orbit is unbound as 21 C is unstable for a neutron emission, and (ii) the 
0^5/2 orbit is bound by at most 2.93 MeV, which is the neutron separation energy of 20 C. 
Without the V s term, the above conditions were barely met only by making V\ much larger 
than the standard strength [l9j. The set A potential in Table U corresponds to this case. 
With the V s term included, we have more freedom to generate different potentials, which 
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TABLE I: Parameters of the n— 20 C potential U . /x = 0.09 fm~ 2 . a and tq are 0.6 fm and 1.3 fm for 
set A, while they are 0.65 fm and 1.25 fm for sets B, C and D. e is the s.p. energy of the n— 20 C 
relative motion. Energy and length are given in units of MeV and fm, respectively. 





V V 1 


Vs 


e(0si/ 2 ) 


£(0p 3 /2) 


e(0pi/ 2 ) 




set A 


33.22 42.10 


0.00 


-19.03 


-9.86 


-4.77 


-1.00 


set B 


43.24 25.63 


9.46 


-19.79 


-14.32 


-11.00 


-2.93 


set C 


41.08 25.63 


7.14 


-19.56 


-12.88 


-9.58 


-1.93 


set D 


38.76 25.63 


4.66 


-19.31 


-11.37 


-8.09 


-0.93 



offer the opportunity of investigating the sensitivity of U on theoretical results. The spin- 
orbit strength V\ is fixed to be the standard value. Three sets of U of this type are listed 
in Table IJ as B, C and D. These potentials are determined by giving different values for the 
0^5/2 single-particle (s.p.) energy: Set B potential gives the deepest energy, while set D the 
shallowest energy. The energies of the lower s.p. orbits turn out to be considerably different. 
It should be noted, however, that our result for 22 C never depends on these energies but on 
their s.p. wave functions as will be seen later. Fortunately, the different potentials chosen 
here give almost the same s.p. wave function for each occupied orbit. All of the potentials 
are set to predict the lsi/ 2 s.p. energy almost zero. It may be probable that the s-wave 
potential strength is further weaker. In that case, the ground state energy of 22 C which we 
will obtain below is to be considered a minimum. 
The ground state of 22 C is described as follows: 

K 

# = $ c $ 2 „, with $ 2n = ^C i $(A i ,A), (3) 

i=i 

where <3> c is the intrinsic wave function of 20 C and the valence neutron part $ 2n is given as 
a combination of correlated Gaussian bases 

$(A, A) = {1- P l2 ) {e-^iy^xMix^Xsih 2)]oo} , (4) 

where Pyi permutes the neutron coordinates and xAx = AuX^+2Ai 2 Xi -x 2 +A 22 x 2 . The 
coordinates a?i=p+|A and x 2 =p — |A are the distance vectors of the neutrons from the 
center of mass (cm.) of 20 C. The angular parts of the two-neutron motion are described using 
yim( r ) — ri Yje m (r) and they are coupled with the spin part xs to the total angular momentum 
zero. The basis function is specified by a set of angular momenta A=(£, S) (L=S), and a 
2x2 symmetric matrix A (A 2 i=A\ 2 ). The two neutrons are explicitly correlated due to 
the term Ay 2 Xi-x 2 , the inclusion of which assures a precise solution in a relatively small 
dimension 

It is vital to take into account the Pauli principle for the valence neutrons in determining 
the energy and corresponding wave function. Though the fulfillment of antisymmetrizing 
the core and valence neutrons is beyond the present model, the Pauli constraint is included 
by imposing that the valence neutrons cannot occupy any s.p. orbits u n ij m of $ c . Here 
Unijm are generated from U, and n£j runs over 0si/ 2 , 0p3/ 2 , 0pi/ 2 , and 0<i5/ 2 . We used 
the stochastic variational method (SVM) to optimize the parameter matrices A. The 
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TABLE II: Properties of 22 C. Length is given in units of fm. 





£(MeV) i?? ms 


BP 


T>m 


M 


) V(p 2 ) 


V(A 2 > 


(Xy X 2 ) Ps=0 


((s 1/2 ) 2 ) 


set A 


-0.413 4.11 


2.44 


3.73 


8.17 


6.19 


10.7 


9.80 0.998 


0.968 


set B 


-0.489 3.96 


2.43 


3.61 


7.54 


5.86 


9.48 


11.9 0.981 


0.915 


set C 


-0.573 3.93 


2.43 


3.58 


7.37 


5.66 


9.44 


9.83 0.990 


0.942 


set D 


-0.388 4.12 


2.44 


3.74 


8.21 


6.29 


10.6 


11.6 0.995 


0.954 



SVM increases the basis dimension one by one by testing a number of candidates which are 
chosen randomly. The basis selection with the SVM is very effective for taking care of the 
short-range repulsion of V\ 2 as well as satisfying the orthogonality constraint. 

The most important channel for the binding of 22 C was found to be A=(0, 0), and other 
channels included were (1,0), (2,0), (1,1), and (2,1). Note, however, that our correlated 
basis functions in practice include higher partial waves as well. Convergent results are ob- 
tained with the basis dimension of 7^^300. The [/-dependence of the solution is moderate 
as shown in Table HT) The result with the ^-independent set A potential is similar to those 
with the other potentials, especially set D potential. This indicates that the present result 
is not very sensitive to the potential provided that it is chosen to satisfy the two conditions. 
The ground state energy is about —390 to —570 keV with respect to the 20 C+n+n threshold, 
which is consistent with the empirical value of —0.423 ± 1.140 MeV [21]. In order to see the 
importance of both spatial and angular correlations of the basis functions, we repeated the 
following calculations. The first was to include only the single channel of A=(0, 0), and then 
the ground state energy turned out to be —0.29 MeV for set B. In the second calculation 
which truncates the basis functions to those with A=(0,0) and A 12 =0 (no correlation cal- 
culation), we obtained the result that the ground state is bound by at most 90 keV. Thus 
the inclusion of the correlated bases is found to gain the energy of about 400 keV. 

The rms neutron, proton and matter radii of 22 C, assuming pointlike nucleons, are listed 
in Table El They are obtained using the corresponding radii of 20 C, 3.23, 2.37 and 2.99 
fm, which are calculated from <3> c . The rms neutron radius is 3.9-4.1 fm. The rms matter 
radius results in about 3.6-3.7 fm, which corresponds to that of a stable nucleus with A ^60. 
Accordingly one may call 22 C a giant halo nucleus. The probability of finding the spin-singlet 
neutrons, Ps=o, shows that the ground state of 22 C almost consists of the 5=0 component. 
Therefore, the non-central potentials have small expectation values and play a minor role in 
binding 22 C: In set B case, the value of (^(tensor + spin-orbit)) is only 7 keV and that of 
(C/i (spin-orbit)) is 57 keV. Thus the binding energy contribution virtually comes from the 
kinetic energy and the central potentials of both U and t> 12 . 

It is interesting to understand how the Borromean system is bound. First of all, we 
note that the non-central forces are found to give negligible contributions. Rewriting the 
kinetic energy as T\ + T p = T\ + T 2 + T IC [22], where Tj is the kinetic energy for the n— 20 C 
relative motion and T rc is the recoil correction term, we decompose the energy contribution 
as follows: 

E = 2(T 1 + U 1 ) + (T TC ) + (v 12 ). (5) 

The decomposition for set B case is 2x(7.185-6.436) -0.118- 1.868 = -0.489 MeV. Except 
for the small contribution of the (T rc ) term, we conclude that the binding of 22 C is obtained 
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by a delicate balance of the two factors: One is that the attraction of vu, though not large, 
keeps the neutrons from separating, and the other is the weak attraction of U which puts 
the neutron in continuum. 

We calculate the probability ((£j) 2 ) of finding the halo neutrons in an (lj) 2 component. 
The ((si/ 2 ) 2 ) value is listed in Table [TTJ Other probabilities are, for set B case, 0.033, 
0.024, 0.009, 0.007, 0.003, 0.003 for £j=d 3/2 , p 3/2 , P1/2, iV/2, 4/2, fs/2, respectively. The 
other potential sets give similar results. We find that the (si/2) 2 component is predominant 
and many other components have small admixtures. Since no bound orbit exists for the 
valence neutron, a realistic shell-model description taking into account these components 
would be hard. On the contrary, the present approach has the advantage that it requires no 
s.p. energies, no matter how high the valence neutrons are excited. 

The halo behavior of 22 C is exhibited through the neutron density, p n (r), which is given 

by 

p n (r) = ($ 2n I p«(|p + r) I $ 2n ) + p h (r), (6) 

where p" (r) stands for the intrinsic neutron density of 20 C, which is calculated from $ c , and 
Ph( r ) is the halo-neutron density with respect to the cm. of 22 C 

2 

Ph(r) = ($2n I S ( X i ~hP~ r )\ $ 2n)- (7) 
i=\ 

The integration of p™ in Eq. © takes care of the fluctuation of the cm. of 20 C around the 
cm. of 22 C. The proton density is given by 

pO(r) = (<S> 2n \(? c (±p + r)\<S> 2n ). (8) 

These densities are displayed in Fig. [TJ The contribution of the halo density to p n exceeds 
that of the core density beyond r=6.2 fm. Note that p\ x (r)/2 is, roughly speaking, the 
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FIG. 2: Comparison of the two-body density distribution of the halo neutrons in C with that of 
the proton and neutron in the deuteron. Set B potential is used. 



squared single halo-neutron wave function. The dip at around r=2.4 fm is due to the 
orthogonality of $ 2n to the 0si/ 2 orbit. 

It is of interest to examine the correlated motion of the two neutrons. A two-neutron sub- 
system with S = is often called a di-neutron when they have spatial extension comparable 
to that of the deuteron. The two-body halo-neutron distribution function, 

p n - n (r) = ($ 2n | 8(X - r) I $ 2n ), (9) 

is compared in Fig. El with the corresponding p — n distribution function of the deuteron 

1 i 

P P -n(r) = - ^ ($d(lM) | 5(r p - r n - r) | $ d (lM)>, (10) 

^ M=-l 

which is calculated using the G3RS potential. It is found that p n _ n (r) has a distribution 
much wider than p p _ n (r). Thus the di-neutron correlation is not prominent in 22 C. The 
value of (T x + v 12 ) is 6.16 MeV (set B), which is to be compared to -2.28 MeV (G3RS) of 
the deuteron. Since p n _ n (r) has a long tail, one may expect that the use of a two-nucleon 
potential with a one-pion exchange tail would give a potential energy different from the 
G3RS potential of a Gaussian tail. To check this point, we estimated the energy difference 
arising when the singlet-even central potential of G3RS is replaced with that of the OPEG 
potential (case 1) [18|, using 



J drp n _ n (r) [ Vl2 (OPEG) - t; 12 (G3RS)] (11) 



and found that the energy gain is only 6 keV. 

Another function of interest is the two-neutron correlation function defined as 

p(x 1 ,x 2 ,9) = ($ 2n | $ 2n ) spin , (12) 

where 6 is the angle between Xi and x 2 and (■ ■ -^pin indicates that the integration is to 
be done over the spin coordinates only. Figure 01 displays the value of 87i 2 x 4 sm9 p(x,x,9). 
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FIG. 3: The two-neutron correlation function p(x,x,6) weighted by 87r 2 2; 4 sin#. The lower panel is 
its contour map. Set B potential is used. 



en 



- 1 1 1 — i — i — r 



1 1 1 1 — i — i — r 

22 p 12 p A 

20 C- 12 C o ■ 

12 c- 12 c 

Ozawa ■ 

Takechi A 

Kox • 




J I I I 1_ 



100 

Incident Energy [A MeV] 



1000 



FIG. 4: Prediction of the reaction cross sections of 22 C and 20 C on a 12 C target. Experimental 



data are taken from Refs. 



One prominent peak appears around £=5.0 fm and 6=17°, which is often attributed to the 
correlation of di- neutron type, but the spatial extension of the two neutrons is too wide to 
be called the di-neutron, as shown in Fig. |21 The peak is followed by a plateau extending 
to larger angles. The valley of the correlation function which appears at around x=2A fm 
reflects the dip observed in the halo-neutron density of Fig. [T] 

The interaction cross section data for the carbon isotopes are available up to 20 C for high 
incident energies With a reaction model proposed in Ref. jlH, we predict the reaction 
cross section cr R of 22 C (and 20 C) using the calculated densities. To make the prediction 
reliable, we modify a nucleon-nucleon (NN) profile function T^tv available in literatures [24| 
so as to reproduce both the elastic scattering cross section and the total cross section of the 
NN collision. Details will be published elsewhere. Figure |U displays <7 R ( 12,20 ' 22 C) on a 
12 C target calculated at several incident energies. A good agreement between theory and 
experiment for ctr( 12 C) confirms the validity of the modification of T^n- The or( 20 C) value 
at the incident energy of 900 A MeV is fairly well reproduced, which indicates that our 
model for $ c is acceptable at least in its prediction for the radius of 20 C. We thus expect 
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that <7r( 22 C), or at least the increase of the cross sections, or( 22 C) — (j R ( 20 C), is predicted 
to good approximation. A measurement of 0r( 22 C) for a wide range of incident energies will 
provide us with valuable information for quantifying the extent to which the halo reaches in 
far distances. We are studying the reaction cross section of 22 C on a proton target as well. 

To conclude, we studied the ground state structure of 22 C in the 20 C+n+n three-body 
model with the orthogonality constraint. The N—1A subshell closure was assumed for 20 C. 
We showed that 22 C is an almost pure ,S-wave two-neutron halo nucleus, and that the non- 
central forces play no active role in binding this fragile system. A measurement of the 
reaction cross section of 22 C+ 12 C is desired to establish the halo structure experimentally. 

We thank A. Kohama for his interest and valuable discussions and M. Fukuda and M. 
Takechi for sending us the reaction cross section data of 12 C+ 12 C. This work was in part 
supported by a Grant for Promotion of Niigata University Research Projects (2005-2007). 
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